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Abstract 

We consider the role played by subcritical bubbles during the electroweak phase 
transition, estimate their average size, amplitude and formation rate taking into 
account the crucial role played by thermalization. We also study the influence of 
subcritical bubbles on the formation of critical bubbles in the thin wall regime 
and show that, contrary to some recent claims, subcritical bubbles do not affect 
the nucleation of critical bubbles in an appreciable way. From this fact we 
conclude that the electroweak baryogenesis scenarios associated with a first order 
electroweak phase transition still remain an attractive possibility. 
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1 Introduction 



Critical bubble nucleation during a first order electroweak phase transition has received 
much attention since the discovery of the possibility for electroweak baryogenesis 
Indeed, one of the basic ingredients for the generation of the baryon asymmetry (apart 
from the requirement of baryon- and CP- violating interactions) is the presence of an 
out-of-equilibrium state |2j which, during the first order electroweak phase transition 
with supercooling, is attained by critical bubbles expanding in the thermal bath of the 
unbroken phase. 

Less attention has been paid to the environment where the critical bubble nucle- 
ation is supposed to occur. Since critical bubbles have a finite size, phase transitions 
are highly local phenomena. Fluctuations with spatial correlations comparable to the 
critical bubble size may be expected to be important for bubble nucleation. Also, 
if thermal fluctuations are too large, any perturbative scheme could break down. In 
such a situation prediction of a first order phase transition becomes suspect, with the 
possibility that the entire scenario of electroweak baryogenesis might be invalidated. 

Although the presence of thermal fluctuations in any hot system is undisputed, 
their role in the dynamics of weakly first order phase transitions is still under debate 

I- 

The idea that statistical fluctuations around equilibrium are spherically symmetric 
and have roughly a correlation volume, where the correlation length is given by the 
inverse temperature dependent mass of the Higgs field, £(T) = m~ 1 (T), was first dis- 
cussed in 0. These fluctuations are referred to as subcritical bubbles. The amplitude 
of the thermal fluctuations was estimated in ref. where it was concluded that they 
are dominant if the Higgs mass tuh is larger than ~ 80 GeV, whence the fraction 
which the asymmetric vacuum occupies in the neighborhood of the critical tempera- 
ture becomes of the order of unity. Therefore it was concluded that critical bubbles 
cannot be generated due to the inhomogeneities of the background field. In ref. |4j], 
however, the continuous disappearance of the subcritical bubbles was not accounted 
for. This can happen in two ways: the subcritical bubbles, being unstable configura- 
tions, tend to shrink; the bubbles are also subject of constant thermal bombardment 
so that they may disappear simply because of thermal noise. The thermalization rate 
of small-amplitude configurations near the critical temperature has been estimated in 
ref. || for the electroweak phase transition, and it was found that, compared with 
typical first order transition times, thermalization is rather fast. 

Kinetics of subcritical bubbles has been investigated by Gelmini and Gleiser |J, 
who found, under a specific assumption about the form of the destruction rate due 
to thermal noise, that for Higgs masses below ~ 55 GeV the approach to equilibrium 
is dominated by shrinking. Unfortunately, for the interesting range of Higgs masses 



dictated by the experimental constraints coming from LEP, m# > 60 GeV J7], their 
analysis is inconclusive since the approximations adopted break down. Nucleation 
induced by the growth of conglomerates of subcritical bubbles has been discussed in 

1- 

Very recently, Shiromizu et al. |J have re-estimated the amplitude of thermal 
fluctuations by calculating the typical size of subcritical bubbles during the electroweak 
phase transition. They claim that for any experimentally allowed values of the Higgs 
mass, the amplitude of thermal fluctuations always exceeds the first reflection point of 
the effective potential. From this analysis they conclude that any standard electroweak 
baryogenesis scenario associated with a first order first transition cannot work. 

Their starting point is the observation that at the microscopic level, the true origin 
of the dominant thermal fluctuation is the perpetual creation and annihilation of 
spherical subcritical bubbles. Thus one should identify the typical size R of the bubbles 
with the size estimated by a statistical ensemble averaging, instead of assuming it to 
be the correlation length. Since the amplitude of the thermal fluctuations sensitively 
depends on the size of the subcritical bubble, a small change in the size results in a 
drastic change in the details of the phase transition. Indeed, modeling any subcritical 
bubbles with a Gaussian shape, the authors of ref. |§ find that the typical size (R) of 
a subcritical bubble is smaller than the correlation length, and that the amplitude of 
thermal fluctuations increases compared with the previous estimates. 

In the present paper we wish to critically reanalyze the results of ref. |J. We 
reconsider the role played by the subcritical bubbles during the electroweak phase 
transition, estimate their size, amplitude and the formation rate. Moreover, we study 
the influence of subcritical bubbles on the formation of a critical bubble by comparing 
the average subcritical energy density associated with the fluctuation to the energetics 
of barrier penetration in the thin wall approximation. Our final conclusion, which 
differs from ref. f9|, is that subcritical bubbles do not affect nucleation of critical 
bubbles in an appreciable way. 

The paper is organized as follows. In Sect. 2 we briefly describe the small super- 
cooling limit and the thin wall approximation. In Sect. 3 we discuss the salient features 
of subcritical fluctuations and present our Ansatz subcritical bubbles. Thermal aver- 
ages are performed in Sect. 4, and the crucial role of thermalization is discussed in 
Sect. 5. Finally, our numerical results and conclusions are given in Sect. 6. 

2 Small supercooling limit and thin wall approxi- 
mation 



First order phase transition and bubble dynamics in the Standard Model have lately 
been studied in much detail, and it has become increasingly clear [[Kj] that for Higgs 



masses considerably heavier than 60 GeV, the electroweak phase transition is only of 
weakly first order. For Higgs mass m# > 100 GeV the calculations, both perturbative 
and lattice ones, confront technical problems^. Therefore, in the paper at hand we use 
a phenomenological Higgs potential for the order parameter <fi suitable for a simple 
description of a first order phase transition: 

V(0) = lm 2 (r)0 2 -iaT0 3 + iA0 4 , (1) 

where we have not determined the parameters perturbatively but fit them, when 
needed, according to a recent two-loop calculation of the gauge invariant effective 
potential flLl|| . Most of the dynamical properties of the electroweak phase transition 
associated with the potential Eq. ([[]), such as the smallness of the latent heat, the 
bubble nucleation rate and the size of critical bubbles, have been discussed in 



For the purposes of the present paper it suffices to recall only some of the results. 

First we need the size of the critical bubble. Assuming that there is only little 
supercooling, as seems to be the case for the electroweak phase transition, the bounce 
action can be written as 

where A(T) = 9Xm 2 (T)/(2a 2 T 2 ). The cosmological transition temperature is deter- 
mined from the relation that the Hubble rate equals the transition rate oc e _5//T , 
yielding S/Tj ~ hi^MpjTj) ~ 150, where Tf is the transition temperature. Thus we 
obtain from Eq. (fj) 

X(T f ) ~ 1 - 0.0442^ = 1-6. (3) 

On the other hand, small supercooling implies that 1 — A = 8 -C 1, i.e. a <C 500A 3 / 2 . 
Solving for A in Eq. (^) yields the transition temperature Tf. One finds 

2o 2 - 

m\T f ) = — \(T f )T 2 . (4) 

The extrema of the potential are given by 

aT 



0±(T) = — (1± V1-8A/9). (5) 



Expanding the potential at the broken minimum 4>+(T) we find 

1 1 a 9 / 2 / x 

e = T f ) = -m 2 (T f )<Pl - -A0 4 + = -0.00218 ^T} + O (S 2 ) . (6) 



1 It is conceivable that for such large Higgs masses the electroweak phase transition is close to a 
second order and does not proceed by critical bubble formation. 



The height of the barrier is situated at ~ 0+/2 with V (0_, T c ) = V max = a 4 T 4 /(144A 3 ), 
where T c is the temperature at which V^(0) = V((f) + ), given by the condition m(T c ) 2 = 
(2 a 2 T 2 /9 A). As T c ~ Tj we may conclude that the thin wall approximation is valid 
if -e/Kiax = 0.314 a 1/2 /A 3/4 < 1, or a < 10A 3/2 . Thus the small supercooling limit 
is clearly satisfied if the thin wall approximation is valid. 

To get the size of the critical bubble we still need the surface tension. One easily 
finds 



2 V2a 3 m3 



* = l dH 2V(T c) = 1 ^Tt. (7) 

We define the critical bubble radius by extremizing the bounce action. The result is 

A 3 / 4 

R C = 13A (8) 
ol I m\± /) 

Therefore R c is much larger than the correlation length £(T/) = l/m(Tj) at the 
transition temperature, as it should. 



3 Subcritical bubbles 



We begin this Section by briefly describing the procedure adopted in ref. || 
following spherical Ansatz was chosen for the subcritical bubble configuration: 



0g 0) = 0+exp 



Rity 



The 



(9) 



where R(t) is the time dependent size of the subcritical bubble and r is the radial 
coordinate. In order to estimate the typical size of a subcritical bubble, one defines 
the canonical momentum pn = dL e g(R, R)/dR, whence an effective Hamiltonian 



H eS (R,p R ) = p R R - L eff 



(10) 



can be constructed. 

The thermal average of R was defined as 



(R) 



J dp R dR R e-P H ° f 
J dp R dR e-^eff 



:i^ 



which, in the neighborhood of the critical temperature, turns out to be much smaller 
than £(T C ) . 

To find the thermal fluctuation amplitude on the scale (R), the following trial 
function was adopted 



(r) = a exp 



(12) 



whence the amplitude in the symmetric vacuum was given by 




J da a 2 e-^ a ^ 
fdae-^ T ) 



(13) 



At T ~ T c the amplitude of a fluctuation y (a 2 ) exceeded the inflection point </>_(T c ) 
and from this fact it was concluded that thermal fluctuations drastically change the 
dynamics of the phase transition from the ordinary first order type with supercooling. 

We now wish to make some critical comments on the approach used in ref. and 
present a remedy to it, which, as we shall argue, is more realistic. 

Let us first make the general observation that it is the actual transition temperature 
Tf rather than the critical temperature T c which is relevant for the study on subcritical 
bubbles. This is true in the sense that if critical bubbles are not important at Tf, they 
most certainly will not be so at T c . As we shall show, it actually turns out that 
subcritical bubbles are not important even at Tf. This justifies, in retrospect, our 
choice T = Tf for performing the calculations. 

In the case of a weak first order phase transition the critical bubble is typically 
well described by a thin wall approximation, where the configuration is by no means 
gaussian, but has a flat 'highland' (with <p determined by the non-zero minimum of the 
potential) and a steep slope down to <p = 0. Therefore it seems natural that also a large 
subcritical bubble should resemble the critical one, i.e. when R increases, the form of 
the subcritical bubble should deform smoothly so that, when R = R c , the bubble is 
a critical one. This behaviour is not reproduced by the Ansatz of Shiromitzu et al. , 
given in Eq. (|9]). Moreover, it has the strange property that in the limit R — > 0, the 
bubble becomes infinitely sharp because the amplitude of the fluctuation is kept fixed, 
while, in general, one would expect that the amplitude of the fluctuations depends 
sensitively on their spatial size. Finally, but no less important, the calculation of 



yj (a 2 ) , Eq. (|T3|), is performed by substituting for the typical size of the subcritical 
bubble the average (R) as calculated from Eq. ([□]), which only is valid for a particular 
value of a, namely a = <p + . In other words, in ref. || the only statistical degree of 
freedom is taken to be the size of subcritical bubbles. 

Motivated by these observations, let us define a subcritical bubble as a functional 
of both the amplitude a and the radius R. For this purpose one has first to study the 
behaviour of the potential as a function of the amplitude. At Tf there is a interval 
£ [a_, a+] where V{<p) < 0. If the amplitude of the bubble is in that interval, 
there exists a critical bubble-solution of the bounce action. This means that we have 
a relation R c = R c (a) which reproduces Eq. (|j) if a — <f> + . Therefore R c (a) serves as 
an upper limit for the bubble radius in that region. 

Moreover, as argued above, the subcritical bubble should be of the thin wall type, 
with the central region having the field value somewhere in the interval £ [a_, a + ]. 




These considerations lead us to define different Ansatze for various regions in the 
(a,R) -plane. When e [a_, a+], we use an Ansatz such that when R — > R c (a), 
the field configuration goes towards the thin wall form. For small R we use a simple 
gaussian configuration. For other values of a we always take a thin-wall like Ansatz. 
Thus we write for e [a_, a + ] and R < R c (a) 

R c — R , i? 



p y p * 



(14) 



where t is the time coordinate^ and 

MR) = e~ r2 / R2 , (15) 
MR) = l/(e m(r - R) + l). (16) 

Such an Ansatz reproduces the requirement that when R — ► i? c , subcritical bubbles 
should resemble critical ones. In practise the statistical averages depend only weakly 
on a because the main contribution to them comes from the region of small a and 
large R. Therefore we assume for simplicity that criticality depends only weakly on a 
and take R c {a) = R c to be a constant whenever possible. 

For <p ^ [a_, a + ] we assume that no gaussian component is present and write 
simply 

4>(t, R) = a(t)M R )- (17) 

However, the statistical averages are excepted to be quite insensitive of the precise 
form of the bubble. 

These Ansatze can be plugged into the action 

S[a, R] = J 'd*x[±(d<f>) 2 -V(<t>)) (18) 

from which the Lagrangian in terms of the dynamical variables a and R can be ex- 
tracted. In the practical calculation we have, whenever possible, approximated <f) t by 
the step function. After that is a simple matter to calculate the effective Hamiltonian 
function H e g of the dynamical variables a and R. 



4 Statistical averages 

Once we have the Hamiltonian, we may calculate the statistical average of a dynamical 
variable of the type F(a, R) simply by 

IF (a R)) = - d?R dPa da dR R " >e ~ PHcB (19) 
/ dpR dp a da dR e~^ eff 



2 Note, however, that we need not to specify the explicit time evolution of a and R when dealing 
with statistical averages. 



However, because the effective Lagrangian is of the form 

£ cS = \(d R)K^ a A-V, (20) 

where K = K(a, R) is a symmetric matrix, after the momentum integration the aver- 
age can be cast into the form 



x v ;/ fdadRVdtiKe-W 

The matrix 

K = aJ Ku M (22) 
V^2i K 22 ) 1 ; 

and the pseudopotential V are given separately for the two regions. For e [a_, a+] 
we obtain 

# 2 K n = A 2 ^ 2 + 2AR A B\ + -i? 5 

3 

i? 2 K 12 = 2aA 2 i? 2 ^ - aAR 3 A\ + ^aR 4 + aR 5 mI(mR) 

+aAR 3 B l 2 - aR 4 B l 2 + 2aAR 3 B\ + aAR 4 mJ 2 {mR) 

R 2 C K 22 = 4a 2 A 2 RAl + a 2 R 5 m 2 I(2mR) +4a 2 AR 3 mJ 4 (mR) 
+a 2 R 3 Aj - 2a 2 R 3 B\ + ^a 2 R 3 - a 2 AR 2 A\ 
-2a 2 R 4 mJ 2 (mR) + Aa 2 AR 2 B{ + 2a 2 R 4 mI{mR) (23) 

V = 2a 2 A 2 iL4 2 + 2a 2 AR 3 mJ 3 {mR) + -a 2 R 5 m 2 I{2mR) 



and 

= 2a 2 A 2 RA 2 , + 2a 2 AR 3 mJ 3 (mR) + - 

+ -m 2 a 2 A 2 i? 3 A 2 + m 2 a 2 A J R 4 5 2 1 + -m 2 a 2 R 5 - -aT-^- 
2 6 3 Racl 



-aT-^A 2 R*B 2 2 - aT-^AR^Bl - laT-^—tf 

Note that in Eq. fl2"3|) the a -dependence of i? c has to be used explicitly because the 
critical behavour is determined from it. For the region where ^ [a_, a + ] the corre- 
sponding functions are given by 

An = i« 3 

K12 = aR 3 mI(mR) 

K 22 = a 2 R 3 m 2 I(2mR) (25) 



A 


= R c {a) - R 






= r duu n e- ku2 = 


r(f) 

2k^ 


K 


Jo 




I(x) 


Jo 


1 2 

X X 2 


J nipt) 


= f 1 duu n e- u2+x{u - 
Jo 





and 

— V = -a 2 R 3 m 2 I(2mR) + -m 2 a 2 R 3 - -aTa 3 R 3 + —\a 4 R 3 . (26) 
4vr 2 v 7 6 9 12 v ' 

A number of shorthand notations have been introduced in the previous equations: 

(27) 
(28) 

(29) 
(31) 

5 Thermalization 

Motivated by the fact that thermal fluctuations can generate configurations with spa- 
tial size comparable to the critical bubble radius, which may affect the dynamics of a 
first order phase transition, the authors of ref. || have estimated the lifetime of fluctu- 
ations of an on-shell Higgs field with zero momentum (po = Tn(T), p = 0). This choice 
reflects the fact that critical bubbles are typically much larger than the interparticle 
distance ~ 1/T in plasma. Writing p = uj — i j/2, one finds that the dispersion 
relation is 

W 2 = |p| 2 + m 2 (T) + i 7 2 , (32) 

where 

ImT( 2 ) , . 

7 = , (33) 

on 

being the two-point function for the Higgs field. 

The imaginary part arises at one loop level, but because of kinematical constraints, 
the two loop contribution is actually dominant in the region of physical couplings. The 
thermalization rate 7 for small amplitude scalar fluctuations and large spatial size, 
R ~ (pp 1 3> 7~\ is estimated [[J to be of the order 7 ~ 1CT 2 T near the critical 
temperature, i.e. much larger than the typical first order transition time. This means 
that all small amplitude fluctuations with size larger than 

i?max = 0(1/7) (34) 



will effectively be absent from the mixture of subcritical bubbles and must not counted 
in the thermal averages. In practise, the limit Eq. (R3) is of the order of few times R c , 
depending on the actual value of 7. Even if it is not precisely known, its inclusion in 



the calculations is important. Without it all statistical averages would be dominated 
by infinite, infinitesimally small fluctuations. Technically this can be seen from the 
Eq. (|2"T|), where the integrals diverge in the limit a — > 0, R — > oo. It is important to 
note that the divergence is not a problem of our Ansatz but merely a more general 
phenomenon, which seems to be related to the general infra-red instability problems 
emerging in the calculations of the effective action. 



6 Results and conclusions 

We have computed the average radius and the amplitude of fluctuations at T = Tf from 
Eq. (|2lD numerically, using a cut-off -R max as discussed in the previous Section. For 
definiteness, we fitted our phenomenological potential Eq. (|I|) to the two loop result for 
the effective potential calculated in |TT[ for the Higgs mass Mh = 70 GeV. This yields 
a ~ 0.048 and A ~ 0.061. One readily verifies that we are indeed safely in the thin 
wall limit. With these parameter values, in units of R c , Tf = 85.70, m 2 (Tf) = 56.78 
and 7 ~ 10~ 2 T f ~ 0.875. 

Choosing, for instance, -R max — 6.6 R c , we find (in units of R c ) 

(a) = 1.16; (R) = 2.93; (a 2 ) = 12.6; (R 2 ) = 12.2 . (35) 

Note that £ = m'^Tf) < (R), so that the average subcritical bubbles are much larger 
than the correlation size, as seems reasonable. Indeed, the correlation length gives the 
order of magnitude of the size of a (thermal) quantum fluctuation, so that, since to 
build up a bubble configuration many quanta are needed, one can expect that (R) 
is larger than £(Tf). Here we remind the reader that the authors of ref. || find the 
opposite result, which is somehow suspect when the above considerations are taken 
into account. 

As far as the production rate for subcritical bubbles is concerned, it is roughly 
given by |§ 

T sc ~ Tf exp (S/Tf) , (36) 



where S, given by Eq. (|18l), is the free energy for the subcritical configuration. 

We have found the average rate for subcritical bubble formation to be much larger 
than the Hubble rate H: 

Y sc = T f (e~ s/T ) ~ 0.49 T f > H . (37) 

Thus we find out that the subcritical bubbling is so fast that effectively the background, 
on which a critical bubble is to formed, is filled by subcritical fluctuations. In principle, 
this could affect the way by which critical bubbles are produced and naively, one 
could think that the extra energy would facilitate barrier penetration or possibly even 



invalidate our thin wall approach. That this is not the case can be ascertained by 
considering the average fluctuation energy about = 0. Using Eq. (|35| ) one obtains 



(V(a)) ~ -m 2 (T f )(a 2 ) = 357.7 < -e ~ 5 x 10 3 , (38) 

where Eq. (|]) and the fact that 0+ ~ 48.05 have been used. Also the average amplitude 
is clearly smaller than the inflection point, \J (a 2 ) <C 0_. Thus, in terms of barrier 
penetration, subcritical fluctuations represent only a minor correction and they have 
no remarkable effect on the thin wall approximation. 

We should point out that the results Eq. fl35|) are rather sensitive to -R max ; for 
instance, if we took -R max — 3.3 R c , we would gotten (a) = 1.94, (R) = 1.53, together 
with an almost twice as large (a 2 ) and an (R 2 ) smaller by a factor of about four. This 
is a general trend: as -R max increases, (a) and (a 2 ) decrease so that the amplitude 
becomes more peaked, while (R) and (R 2 ) grow. To explain analytically, for instance, 
the scaling of {R) with the changing of -R max one can consider the limit R c — > oo in 
Eq. (PIP which is reasonable in the vicinity of T c . In such a case, it is easy to show 
that 

J Rmax dR R 1 ' 2 

{ } * f R ™*dRR-V 2 ~ max ' (39) 

whereas a similar calculation proves that (R 2 ) oc -R max . It is then clear that why, 
increasing i? max by a factor two, (R) and (R) 2 increase by a factor two and four, 
respectively. Analogous considerations explain the decreasing of (a) and (a 2 ) with 
increasing R max . 

Note, however, that, for reasonable values of the cut-off R max , the subcritical nu- 
cleation rate does not change, nor the conclusion that (V(a)) <C — e. 

To conclude, we have find that subcritical bubbles are not, at least for Higgs masses 
less than about 100 GeV, important during the onset of the electroweak phase transi- 
tion. We have argued that critical bubble formation proceeds essentially unchanged, 
and here we disagree with the result of Shiromizu et al. |§ . The main reason for this 
discrepancy is that in ref. |§ the only statistical degree of freedom describing the prop- 
erties of subcritical bubbles is taken to be their size R, keeping their amplitude fixed at 
the value + . This fact led the authors of ref. |§ to conclude that (R) <C £(T C ) which 
appears not reasonable. Indeed, to build up a classical configuration, many thermal 
quanta (whose typical size is the correlation length) are needed and one would then 
expect that (R) ^> £(T C ). To recover such a behaviour, we have treated both the 
amplitude and the critical size of subcritical bubbles as statistical degrees of freedom. 
We have then estimated their average size, amplitude and formation rate taking into 
account the crucial role played by thermalization, showing that subcritical bubbles do 
not affect the nucleation of critical bubbles in an appreciable way. From this fact we 



conclude that the electroweak baryogenesis scenarios associated with a first order elec- 
troweak phase transition remain a viable possibility to explain the primordial baryon 
asymmetry in the Universe. 
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